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ABSTRACT 

We examine the prospects for raeasuring the dark energy equation of state parameter w within the con- 
text of any uncertain redshift evolution of galaxy cluster structure (building on Majumdar & Mohr 2003) 
and show that including the redshift averaged cluster power spectrum, {Pcl)^ and direct mass measure- 
ments of 100 clusters helps tremendously in reducing cosmological parameter uncertainties. Specifically, 
we show that when combining the redshift distribution and the power spectrum information for a par- 
ticular X-ray survey (DUET) and two SZE surveys (SPT & Planck)^ the constraints on the dark energy 
equation of state w can be improved by roughly a factor of 4. Because surveys designed to study the 
redshift _distribution of clusters will have all the information necessary to construct Pcu the benefit of 
adding Pd in reducing uncertainties comes at no additional observational cost. Combining detailed mass 
studies of 100 clusters with the redshift distribution improves the parameter uncertainties by a factor of 
3-5. The data required for these detailed mass measurements- assumed to have la uncertainties of 30%- 
are accumulating in the the XMM-Newton and Chandra archives. The best constraints are obtained 
when one combines both the power spectrum constraints and mass measurements with the cluster red- 
shift distribution; when using the survey to extract the parameters and evolution of the mass-observable 
relations, we estimate the uncertainties on w of ~4% to 6%. These parameter constraints are obtained 
from self-calibrating cluster surveys alone. In combination with CMB or distance measurements that 
have different parameter degeneracies, cluster studies of dark energy will provide enhanced constraints 
and allow for cross-checks of systematics. 

Subject headings: cosmic microwave background — galaxies: clusters — cosmology: theory 



1. INTRODUCTION 

Galaxy clusters have been used extensively to determine 
the cosmological matter density parameter and the ampli- 
tude of density fluctuations. The number density of mas- 
sive clusters is exponentially sensitive to the amplitude of 
initial Gaussian fluctuations [whose normalization we de- 
scribe using as, the rms fluctuations of overdensity within 
spheres of S/i^^ Mpc radius] (Eke et al. 1996; Henry 1997; 
Bahcall & Fan 1998; Viana & Liddle 1999; Reiprich & 
Bohringer 2002; Lin & Mohr 2003). In principle, clusters 
contain information about many aspects of the evolution 
of the universe. Wang & Steinhardt (1998) argued that 
a measurement of the changes of cluster number density 
or abundance with redshift would provide constraints on 
the dark energy equation of state parameter w = p/p. 
Haiman et al. (2001) showed that with future large sur- 
veys it should be possible to obtain precise measurements 
of the amount f2 e and nature w of the dark energy. Since 
then, a series of papers has focused on the prospects of 
high yield cluster surveys as probes of dark energy (Holder 
et al. 2001; Weller et al. 2001; Levine et al. 2002; Hu & 
Kravtsov 2003; Majumdar & Mohr 2003; Hu 2003). These 
papers also show that constraints from cluster surveys will 
be complimentary to those from cosmic microwave back- 
ground (CMB) anisotropy and SNe la distance measure- 
ments. 

The potential of using cluster surveys can only be re- 
alized fully if we not only know the cluster redshifts but 



also have some estimate of the cluster masses (Haiman 
et al. 2001). However, this total mass has to be inferred 
indirectly through use of cluster scaling relations and ob- 
servables like the X-ray flux, SZE flux, galaxy light and 
weak lensing shear. Recently, Levine et al. (2002) exam- 
ined an X-ray cluster survey, showing that a sufficiently 
large survey allows one to measure cosmological param- 
eters and constrain the cluster mass-observable relation 
simultaneously if one assumes perfect knowledge of the 
redshift evolution of galaxy cluster structure. If one al- 
lows for an imperfect understanding of galaxy cluster evo- 
lution, then the constraints on w are seriously weakened 
(Majumdar & Mohr 2003, hereafter MM03). To over- 
come this important limitation, MM03 examined the im- 
pact of follow-up of a fraction of the survey cluster sample 
and showed that even a 1% follow-up helps in improving 
the w constraints by a factor 2 to 3. These calculations 
underscore the importance of incorporating information 
from multiple observables into future cluster surveys, and 
they demonstrate that cluster surveys are essentially self- 
calibrating- containing enough information to solve for the 
mass-observable relation at every redshift and constrain 
cosmological parameters. As another solution to overcome 
the evolution-cosmology degeneracy, Hu (2003) has shown 
that there is potential for internal calibration by using the 
shape of the mass function of clusters in redshift slices. In 
the case of tight mass-observable relations and accurate 
redshifts, this approach essentially enables one to solve for 
the mass observable relation independently in each redshift 
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slice of the survey. Additional calculations are needed to 
examine how well this very promising approach will work 
in future cluster surveys. 

Here we look at another way to improve constraints on 
the equation of state of dark energy. We show how one may 
use the redshift averaged cluster power spectrum Pci{k), 
derived from the correlated positions of galaxy clusters 
in the survey, to tighten the constraints on cosmological 
parameters. Because the sky position and the redshift (al- 
ready needed to measure the redshift distribution) for each 
cluster gives us enough information to construct Pd, one 
essentially gets the cluster power spectrum without further 
observational effort. The amplitude and shape of the clus- 
ter power spectrum depends on the underlying dark matter 
power spectrum and the bias of the tracers, which is re- 
lated to the mass of the halo (Mo & White 1996; Sheth & 
Tormen 1999) in the same way as the mass function of the 
halos. Thus, a measure of biased cluster power spectrum 
essentially provides an indirect measurement of the mass 
of the cluster tracers . This additional information is es- 
sential in constraining the evolution of the mass-observable 
relation, which is key to obtaining tight cosmological con- 
straints from the cluster redshift distribution. In addition, 
the shape of the cluster power spectrum provides an addi- 
tional constraint on the matter density because of its effect 
on the transfer function. 

Moreover, an advantage of Pd over any spatial correla- 
tion, ^(r), studies is in the respective imcertainty analy- 
sis. For Pel, one can assume, for an initial Gaussian field, 
the statistical fluctuation to be mainly uncorrelated (Feld- 
man et al. 1994). However, Pd is similar to ^(r) in the 
sense that both exhibit different parameter degeneracies 
than those from using dN/ dz and hence combining power 
spectrum constraints with the cluster redshift distribution 
helps in tightening parameter constraints. The aim of this 
paper is to examine the complementarity of the redshift 
distribution and cluster power spectrum that can be ob- 
tained in any high yield cluster surveys. We show that by 
combining these two analyses, dramatic improvement on 
w-constraints can be realized even in the face of uncertain 
evolution of the cluster mass observable relation. 

The paper is arranged in the following way. In §2 we 
describe the basic formalism and then in §3 we describe 
three representative surveys and our method of estimat- 
ing parameter constraints. In §4 we present our results 
and finally conclusions are drawn in §5. 

2. PRELIMINARIES 

2.1. Cluster Redshift Distribution 

To begin, we estimate the redshift distribution of de- 
tectable clusters within a survey solid angle AQm, 



dN dV ^ ^ 



dz 



dzdfl 



m,z)'-^dM (1) 



where dV/dzdfl is the comoving volume clement, 
{dn/dM)dM is the comoving density of clusters of mass 
M, and /(M, z) is the redshift dependent cluster selection 
function for the survey. In this analysis we take /(M, z) 
to be a step function at some limiting mass Mum, which 
corresponds to the mass of a cluster that lies at the sur- 
vey detection threshold. In a real survey the selection 
function will depart from a step function, and this has 



been examined elsewhere (Holder et al. 2000; Lin & Mohr 
2003). We use the cluster mass function dn/dM deter- 
mined from structure formation simulations (Jenkins et al. 
2001). Sample dn/dz plots can be found in MM03. 

For an X-Ray survey, the mass-observable relation takes 
the form of bolometric flux - mass relations and we use 

f,{z)4Trdl = A^M^o^E^z) (1 + zf^ (2) 

where fx is the observed flux in units of erg s~^cm~^, 
c?L is in units of Mpc, M200 in units of IO^^Mq is the 
mass enclosed within a radius r2oo having a overdensity of 
200 with respect to critical and H{z) = HqE{z). where 



3(l+to) 
E 



We con- 



E'^{z) = viuii + zf + f7fe(i + zf + n 

vert M200 to M(z), the halo mass appropriate for our 
mass function at redshift z using a halo model (Navarro 
et al. 1997, hereafter, NFW). Following MM03, the pos- 
sibility of non-standard evolution of the mass-observable 
relation is included through the parameter ■ In our fidu- 
cial model we take 7^ = to be consistent with the ob- 
served weak evolution in the luminosity-temperature re- 
lation (Vikhlinin et al. 2002), and we choose (3x = 1.807 
and log Ax = -3.926 (Reiprich & Bohringer 2002). 

For an SZE survey, the corresponding SZE flux-mass 
relation is given by 

/,,(z, v)d\ = f{iy)ficMA,,M^^SE^/^{z) (1 + z)^" (3) 

where f{v) is the frequency dependence of the SZE dis- 
tortion, fsz is the observed flux in mJy, M200 is in units 
of IO^^Mq, ficM = 0.12 (e.g Mohr et al. 1999) and dA is 
in units of Mpc (see Diego et al. 2001). For our flducial 
model, we use the mass - temperature relation for clus- 
ters having T > 3KeV (Finoguenov et al. 2001), to get 
log Asz = 16.47, f3sz = 1.68 and jsz = to model stan- 
dard structure evolution. Note that in determining the 
estimated uncertainties on cosmological parameters, we al- 
low the normalization of these mass-observable relations 
to be free to vary. Finally, we impose a minimum cluster 
mass of 10^^/i~^ Mq on all surveys. 

2.2. The Cluster Power Spectrum 

The power spectrum of future high yield surveys pro- 
vides a different probe of the background cosmological 
model. Taken alone, the constraints from cluster power 
spectrum are far weaker than those from dN/dz alone. 
However, when combined with other measures of cosmo- 
logical parameters, the power spectrum studies can give 
very powerful constraints. While in most cases, people 
have studied the power spectrum for galaxy surveys, the 
cluster power spectrum holds similar promise because clus- 
ters arc highly biased when compared to galaxies, making 
it possible to obtain similar statistical uncertainties with 
far smaller samples. Another important issue that con- 
trasts the cluster power spectrum, Pd{k), from that of 
galaxies is the fact that the bias, b{M, z), for clusters can 
be determined from large scale N-body simulations and 
theoretical calculations (Mo & White 1996; Sheth & Tor- 
men 1999). In addition, cluster masses are related to sim- 
ple observables, making it possible to directly connect the 
bias of the cluster power spectrum to these same observ- 
ables. The same physics of gravitational collapse of dark 
matter dominated halos underlies the theoretical and nu- 
merical predictions of the bias-mass relation and the mass 
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function of collapsed haloes (Mo & White 1996; Sheth & 
Tormen 1999; Jenkins et al. 2001; Sheth et al. 2001). 




Fig. 1. — The upper panel shows redshift averaged cluster power 
spectrum, P^.;, for the SPT survey, for the three redshift bins, 
z = - 0.42 (short-dashed Une), z = 0.42 - 0.72 (long-dashed) 
and 2 = 0.72 — 3 (dotted line), each bin having ~ 9700 clusters. 
The solid line shows the results for one single redshift bin from 
2 = — 3. The lower panel shows the fractional uncertainties in 
estimating the power spectrum, up ^(k)/ P^i{k), using the method 
detailed in Feldman et al. (1994). 

To calculate the parameter sensitivity of Pci{k)we fol- 
low the method described in a number of papers dealing 
with galaxy power spectrum (Tegmark 1997; Eisenstein & 
Hu 1999), keeping in mind that clusters are biased with 
respect to the underlying matter distribution. We use the 
bias model in Sheth & Tormen (1999). In the limit where 
the survey volume is much larger than the scale of any 
features in the power spectrum (Hamilton 1997), the un- 
certainties can be well estimated with the method given by 
Feldman et al. (1994, hereafter FKP). The cluster power 
spectrum is constructed as Pci{k,z) = h\f j{z)P{k, z) 
where P{k,z) is the matter power spectrum and 6e// at 
a particular redshift is the cluster mass function averaged 
linear bias, b{M, z) (Suto et al. 2000; Moscardini et al. 
2000). The effective bias is, then, given by 



JdM^b{M,z) 



(4) 



where dn/dM is the cluster mass function obtained from 
N-body simulations (Sheth & Tormen 1999; Jenkins et al. 
2001). This is valid for a small redshift interval Sz at a 
particular z. It is useful to measure the power spectrum 
over a wide redshift range. In this case, the light cone av- 
eraged cluster power spectrum, Pd, in the redshift range 
of interest between Zmin and z^ax is obtained as 

Pciik) = ^'-'\ '':\ (5) 



where the comoving number density of clusters in the flux 
limited surveys is computed by integrating Eqn 1. Thus, 



Pel not only depends on cosmology but also on the spe- 
cific survey that determines the redshift dependence of ob- 
served clusters. 

We assess the sensitivity of Pd to different parame- 
ters through the Fisher matrix formalism detailed in the 
next section. Note that the evolution of b{M, z) follows 
the growth of structures in hierarchical models of struc- 
ture formation (Mo & White 1996; Sheth & Tormen 1999; 
Sheth et al. 2001). Also note that the amplitude of the 
power spectrum evaluated on the light cone for a given 
flux limit increases with z; the increasing bias at higher 
redshift offsets the reduction in amplitude of the under- 
lying power spectrum of dark matter fluctuations. For 
all the surveys, the power spectrum is calculated in the 
fc-range 0.0014ft.-iMpc"i to 1.4/i~iMpc"^ The uncer- 
tainty in the estimation of Pd increases rapidly as we de- 
crease A;, and effectively cuts off the contribution to the 
Fisher matrix for wavenumbers less than kmin- At higher 
fc, non-linear effects become important and increase Pd- 
Although, the power spectrum in the non-linear regime 
remains a well defined quantity, it is unwise to to use the 
Fisher matrix approach (see below) to estimate constraints 
on parameters as the density field becomes nonlinear and 
non-Gaussian; in this regime Equation 7 could give in- 
correct estimates of the parameter constraints. One needs 
to use Monte-Carlo methods to accurately study processes 
having non-Gaussianities. We ask the reader to keep these 
considerations in mind when assessing our results. Finally, 
we neglect redshift space distortion and non-linear contri- 
butions to Pd , which produce a small change in the power 
spectrum amplitude (Suto et al. 2000). 

2.3. Follow-up of Large Solid Angle Surveys 

Finally, MM03 introduced the idea of follow-up of a frac- 
tion of the cluster sample to calibrate the mass-observable 
relation and highlighted the resulting improvement in pa- 
rameter constraints. FoUowup requires detailed X-ray, 
SZE or galaxy spectroscopic observations that allow one to 
measure the mass-like quantity Mf{9) = M{0)/dA, which 
we will refer to as the follow-up mass, where the M{d) can 
be calculated using a NEW profile for the matter distri- 
bution. We assume a concentration c — 5. At fixed red- 
shift, this follow-up produces an Mj-fx of Mf-fsz relation 
which would provide direct constraints on the structural 
evolution of the clusters. The parameter sensitivity of 
these scaling relation observations can exhibit quite differ- 
ent degeneracies than for the cluster redshift distribution, 
making the scaling relations and the cluster redshift dis- 
tribution complementary (see MM03). We emphasize that 
a large number of clusters are now being studied in great 
detail with Chandra, XMM and deep optical weak lensing 
data. This growing ensemble of well studied clusters will 
be available for mass calibration in all future surveys. We 
model this additional information by placing 100 clusters 
uniformly over a range of redshifts from 0.3 < z < 1.2 and 
masses from lO^'^h'^MQ to lO^^h'^MQ. 

3. COSMOLOGICAL SENSITIVITY OF A SURVEY 

3.1. Fisher Matrix Technique 

We employ the Fisher matrix technique to probe the rel- 
ative sensitivities of cluster surveys to different cosmolog- 
ical and cluster structural parameters. The Fisher matrix 
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information for a data set (see Tegmark ct al. 1997; Eisen- 
stein et al. 1998) is defined as Fij = — < >, where 

C is the likehhood for an observable (proportional to ^ 

for the redshift distribution, P^i for the power spectrum 
and Mf for the follow-up) and pi describes our parameter 
set. The inverse describes the best attainable co- 
variance matrix [Cij] for measurement of the parameters 
considered. The diagonal terms in [Cij] then give the un- 
certainties on each of our parameters. To estimate the full 
potential for cluster surveys, we calculate the final Fisher 
matrix is the sum of the Fisher matrices for the cluster red- 
shift distribution (Ffj), the cluster power spectrum {F^j), 

the follow up observations {f/j) and the priors. 

We construct the survey Fisher matrix Ffj following 
Holder et al. (2001) as 

where wc sum over n redshift bins of size Az — 0.01 to 
Zmax = 3.0 and Nn represents the number of surveyed 
clusters in each redshift bin n. The Fisher matrix for the 
power spectrum is given by 

dPcijk) dPajk) Akm 

™ dpi dpj a'^FKP 

where wc sum over m logarithmic bins in /c-space between 
kmin and kmax (scc, for example Tegmark (1997)). The 
variance in each /c-bin is calculated using the formulation 
given in Feldman et al. (1994) with optimal weighting for 
each survey. 

The Fisher matrix for the follow-up is constructed as 



F^- 



dMf dMf 1 

dpi dp, a 



(8) 



Mf ^ 



where Mf is a function of halo mass M and angular radius 
^ at a particular redshift z. As noted above, we model the 
available foUowup information from deep X-ray studies of 
clusters by summing over 100 clusters that uniformly span 
a range of redshifts (0.3 < z < 1.2) and a range of masses 
(IO^/i-^Mq to IO^^/i-^Mq). In this way we attempt to 
model the information provided by all the detailed cluster 
analyses that acciimulatc in the literature. Any foUowup 
of clusters in a particular survey would result in additional 
constraints on cluster structural evolution. In this analysis 
wc take (Tm, = O.SAfj to be the characteristic uncertainty 
in the follow-up mass measurements. For more details on 
cluster follow-up specifications, please refer to MM03. 

3.2. Future Surveys, Fiducial Cosmology and External 

Priors 

For our study, we adopt three representative surveys, 
two in SZE and one in X-Ray. For the X-Ray survey 
wc examine the 10** deg^ flux limited X-ray survey pro- 
posed as part of the DUET mission to the NASA Medium- 
class Explorer Program. For the SZE surveys, we con- 
sider the 4,000 deg^ SZE survey to be carried out with 
an 8m telescope being constructed for the South Pole 
(SPT) and the Planck full sky SZE cluster survey. We 
model the DUET X-ray survey as having a bolometric 
flux limit of fx > 1.25 x 10~^^ erg/s/cm^ (corresponding 



to fx > 5x 10^^^ erg/s/cm^ in the 0.5:2 keV band). For 
our fiducial cosmological model, this survey yields ^22,000 
detected clusters, consistent with the known X-ray log N- 
logS relation for clusters (e.g. Gioia et al. 2001). We 
model the SPT SZE survey as a flux limited survey with 
fsz > 5 mJy at 150 GHz. Within our flducial cosmologi- 
cal model this survey would yield ~ 29, 000 clusters with 
measured fluxes. And finally for the Planck survey, we 
take a flux limit of fsz > 50 mJy at 353 GHz (at which 
the angular resolution of Planck would reach 5 arc min) 
(Diego et al. 2002) which for 70% of the sky yields ~2 1,000 
clusters. 

The fiducial cosmological parameters (and priors) of our 
model are adopted from the first year WMAP result (Ben- 
nett ct al. 2003) and are h = 0.71, f^M = 0.27, n = 0.93, 
= 0.044, 0-8 = 0.84 and w = -1. From WMAP results 
we know that the curvature is very close to zero, and hence 
wc assume a flat universe, taking Qtot = + = 1- 
We also neglect any variation in w with redshift. The pri- 
ors used for our calculations are ah = 0.035, an^ = 0.004 
and cr„ = 0.03. We do not use any priors on Qm , cr^ 
or w for the main results of this paper as summarized in 
Table 1. However, to show the effect of strong CMB pri- 
ors on , cg and w (shown within parentheses in Table 
1), we use a^M = 0.04, cjcrg = 0.04 and = 0.11. Note 
that this provides a crude estimate of the complementarity 
with CMB constraints, because we have not used the pa- 
rameter covariance appropriate for CMB power spectrum 
analyses. 

For power spectrum calculations, each survey is divided 

into three redshift bins such that each bin has an ap- 
proximately equal number of clusters. For example, for 
the SPT survey which has the most clusters at high red- 
shift, the redshift binnings are ^ = 0.0 - 0.42, 0.42 - 0.72 
and 0^72 — 3.0. For the fiducial cosmology and SPT sur- 
vey. Pel and the the fractional imccrtainty, crpKP / P{k), 
for the three redshift bins are shown in Figure 1 using 
short-dashed, long-dashed and dotted lines respectively. 
The solid line is for the case of a single redshift bin from 
2; = 0.0 -3.0. 

4. COSMOLOGICAL PARAMETER CONSTRAINTS 

Our main results are summarized in Table 1 and high- 
lighted in the following three figures, one for each survey. 
The table contains the l-cr parameter uncertainty for each 
survey for f^M , and w and addresses the main focus of 
this paper, which is to examine the degree to which one can 
tighten cosmological from the redshift distribution dn/dz 
by adding the cluster power spectrum Pd and detailed fol- 
lowup of 100 clusters. With the exception of Majumdar 
& Mohr (2003) and Hu (2003), previous studies of cluster 
surveys have implicitly assumed full knowledge of cluster 
structure evolution and only addressed the the sensitivity 
of cluster surveys to cosmological parameters. Constraints 
from such considerations are grouped under the heading 
'Only Cosmology' in Table 1. However, if one does not 
assume perfect knowledge of cluster structure and lets the 
survey solve for both cosmology and cluster scaling pa- 
rameters, then the cosmological constraints are weakened. 
This is illustrated in the next column labeled 'Std Evol.'. 
This approach, however, still assumes that there is no un- 
expected evolution in scaling relation with redshift and 
ignores any uncertainty due to it. A cautious approach 
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Table 1 

Estimated Parameter Constraints 



Survey /Parameter 


Only Cosmology 


Std Evol. 


Self-Cal 


Self-Cal 


Self-Cal 


Self-Cal 




dfi/dz 


dfi/dz 


dn/dz 


dn/dz-h- P^i 


dxi /d2:+Fup 


dxi /d2:+Both 


.'JPT ,1ZE f^lir'iiP'ii 

kJJ. -L I^ZJUI kJ Uil UCU 
















0.0084 


0.0190 


0.0237 


0.0115 


0.0181 


n 01 08 ("0 0096"! 


Acts 


0.0105 


0.0236 


0.0441 


0.0196 


0.0218 


0.0125 (0.0111) 


Aw 


0.0371 


0.0663 


0.1745 


0.1543 


0.0602 


0.0376 (0.0335) 


Planck SZE Survey 














AHm 


0.0078 


0.0575 


0.0820 


0.0188 


0.0344 


0.0062 (0.0059) 


Acts 


0.0123 


0.1097 


0.1747 


0.0393 


0.0629 


0.0098 (0.0094) 


Aw 


0.0397 


0.1977 


0.3930 


0.1049 


0.1149 


0.0404 (0.0390) 


DUET X-Ray Survey 














AQm 


0.0099 


0.0290 


0.0372 


0.0114 


0.0142 


0.0096 (0.0083) 


Aas 


0.0131 


0.0503 


0.0846 


0.0127 


0.0182 


0.0114 (0.0100) 


Aw 


0.0505 


0.0806 


0.3540 


0.2065 


0.0765 


0.0625 (0.0527) 



would be to try to solve for any evolution from the sur- 
vey itself (see Eqn 2 and 3). The estimated uncertainties 
for this self-calibrating case are shown in column labeled 
'Self-Cal'. The next column shows constraiiits when one 
combines information from both dn/ dz and Pd , both cal- 
culated for the self-calibration case. The second last col- 
umn shows the constraints if one does a mass follow-up of 
100 clusters from the survey sample and the final column 
shows the cosmological constraints when one uses informa- 
tion from dn/dz, Pd and also follow-up. The constraints 
are using cluster surveys only. However, in the last col- 
umn, within parentheses, we show an estimate of how the 
constraints change if one includes strong priors of CIm , 
and w. 

From the table, one is able to draw two general conclu- 
sions: 

• Assuming complete knowledge of cluster structure 
can give overly optimistic constraints on cosmo- 
logical parameters. A complete lack of knowledge 
of cluster structure can drastically degrade our 
ability to constrain such parameters (especially w) 
from future cluster surveys. This is illustrated in 
columns 2-4. 

• The last three columns show that one can greatly 
improve the constraints if one uses information 
from power spectrum studies or mass follow-up of 
100 clusters. Using all available information from 
dn/dz, Pd and follow-up can deliver uncertainties 
comparable to the case of using the redshift 
distribution alone and assuming perfect knowledge 
of cluster structure and evolution. This, indeed, is 
very encouraging. 

4.1. Self- calibration and degradation of constraints 

For the case with perfect knowledge of cluster structure, 
all the three surveys compare favourably, yielding 1-cr un- 
certainty on w of 0.0371, 0.0397 and 0.0505 for the SPT, 
Planck and DUET surveys respectively. The marginally 
better uncertainty on w from SPT over Planck is due to 
the higher sensitivity of the SPT survey, which delivers 
higher redshift clusters that are more sensitive to change in 
w (for example, see Figure 4 in MM03). Planck, which es- 
sentially detects only high mass clusters over the entire ex- 
tragalactic sky, constrains Om most effectively. Note, that 



when we have no uncertainties entering from any imprecise 
knowledge of cluster structure, SPT having the maximum 
number of clusters and the deepest redshift range best con- 
strains as- However, with the addition of extra informa- 
tion from the cluster power spectrum and mass follow-up, 
the shallower but largest solid angle Planck survey pro- 
vides the the best constraints on cts- The joint 1-a con- 
straints for this case on w and are shown in Figure 2 
and on ug and flM in Figure 3 using dotted lines. 

When one solves for both the normalization and the 
slope of cluster scaling relations, along with cosmology, 
from cluster surveys, the uncertainties degrade by factors 
of upto 2-5 depending on the survey and the parameter 
concerned. The fact that one still obtains uncertainties on 
w of 0.0663 and 0.0806, for SPT and DUET respectively 
shows that it is possible to constrain both cosmology and 
cluster physics reasonable well from future cluster surveys 
(see also Lcvinc et al. 2002). The weakening of param- 
eter constraints is most significant for Planck for which 
the peak of dN/dz falls at z ^ 0.14 which is much be- 
low the modal redshift in SPT or DUET. This is due 
to fact that lower sensitivity surveys that are picking out 
only the few, most massive; clusters can less effectively con- 
strain the cluster mass-observable relation normalization 
and slope (i.e its more costly to solve for cluster structure) 
and these larger uncertainties lead to weakened cosmolog- 
ical parameter constraints (for example, Planck can con- 
strain the slope of the mass-observable relation to ~ 27% 
compared to ~ 17% and 9% from DUET and SPT). 

The crux of the problem arises when one takes into ac- 
coimt the possibility of non-standard evolution. When one 
requires that the survey also solve for the evolution of the 
mass-observable relation, the constraints on w weaken fur- 
ther by a factor of ~ 2 — 2.5 for the two SZE surveys 
(i.e ~ 0.1745 for SPT and ~ 0.3930 for Planck) and 
a factor ^ 4.5 for DUET. Because this imcertain evo- 
lution affects the mass estimates for high redshift clus- 
ters, it affects the constraint from DUET more than those 
from Planck, even though both have similar total num- 
ber of clusters. Constraints on Q.m and erg also weaken 
significantly. The reason for this weakening of constraints 
is the fact that constraints on the evolution parameters 
{ix/lsz) are very weak and this large uncertainty in the 
evolution of the mass-observable relation makes it harder 
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to know the masses of high redshift clusters; this drives the 
weakened sensitivity to other cosmological parameters (for 
additional discussion, see MM03). The interplay between 
redshift depth and total number of clusters is complex, 
though, for SPT, though deeper than DUET, has slightly 
less degradation in errors due to it having ~ 30% more 
clusters than DUET. In Figures 2 and 3, we have shown 
constraints on w — ^Im and erg — Qm for 'Self-Cal' using 
long-dashed lines. 
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Fig. 2. — Constraints on w and Qjvf for an SZE SPT survey (top), 
SZE Planck survey (middle) and an X-ray DUET survey (bottom). 
Contours denote joint 1— cr constraints in five scenarios: constraints 
from dN/dz for 'Only Cosmology' case (i-dotted) ; constraints from 
dN/dz for_'Non-Std Evol' case (ii-long-dashed); constraints from 
dN/dz + P^i, (iii-dot-dashed); dN/dz + 100 cluster follow-up (iv— 
short-dashed) and the finally constraints from dN/dz+P^i + follow- 
up (v— solid). Note, that for the Planck survey, for the 'Non-Std 
Evol' scenario (i.e long-dashed), the constraint ellipse is too large 
and lies outside the range for the two axes. A flat universe is as- 
sumed for all the cases. The different terms are described in the 
text. 

4.2. Including the Cluster Power Spectrum 

The addition of constraints from Pd from the survey 
sample to that from dn / dz provides us with a rather neat 
and effective option to dramatically improve the cosmolog- 
ical constraints. The fact that parameter degeneracies on 
cosmological parameters from the redshift distribution and 
the spatial two point correlation function ^ (r) are quite dif- 
ferent has been demonstrated for future surveys (Refregier 
et al. 2002). Because P{k) and ^(r) are Fourier transform 
pairs, this also holds for dn/dz and Pd as well. Includ- 



ing the cluster power spectrum helps in two ways. First, 
because the amplitude of the cluster power spectrum de- 
pends on the cluster bias, a measure of the power spectrum 
essentially provides a handle on the masses of the tracers 
used to measure the power spectrum. This information 
provides an indirect measurement of of masses of the clus- 
ters for any particular survey and mass-observable rela- 
tion. Secondly, the shape of the power spectrum reflects 
the transfer function for density perturbations, which is 
very sensitive to the epoch of matter-radiation equality; 
thus, the cluster power spectrum places strong constraints 
on cosmological parameters, mainly ^m- Because of the 
complex parameter degeneracies, improved constraints on 
VLm are reflected by tighter constraints on other parame- 
ters. 

We emphasize that any survey that delivers the red- 
shift distribution naturally gives the spatial position and 
the redshift of the tracers; hence the cluster power spec- 
trum comes without any extra observational effort. Al- 
though serendipitous cluster surveys may be adequate for 
measuring the cluster redshift distribution, it will require 
contiguous coverage to keep survey window functions suf- 
ficiently simple to make robust estimates of the cluster 
power spectrum. 

It is clear from Table 1 that adding the Pd constraint re- 
duces uncertainties on w for the two larger area surveys by 
factor of ~ 2 — 4 to make the uncertainty closer to ones ob- 
tained from the redshift distribution in the 'Only Cosmol- 
ogy' scenario. Using the cluster power spectrum tightens 
the constraint on w most for Planck, followed by DUET 
and SPT. This is expected because of the large physical 
scales probed by the all sky Planck survey. This is more 
evident when one compares DUET with Planck, both of 
which have similar total number of clusters. Thus, when 
one combines Pd information to that from dN/dz, one ob- 
tains 1-cr constraint of ~ 10% on w from Planck which is 
a factor of ~4 improvement on the ~ 39% constraint on 
w obtained in the 'Self-Cal' scenario. Notice that addition 
of the power spectrum also very effectively tightens the 
constraint on JImj especially for Planck, which is reduced 
by a factor of - 4.5 from crn„ = 0.082 to 0.0188. For SPT 
and DUET, the Q.m constraint also improves by factor of 
- 2-3. 

A point to note is that the amplitude of the power spec- 
trum, which depends on the bias of the tracers and the 
underlying amplitude of the dark matter power spectrum- 
fjg" is a key component which helps in tightening the con- 
straints. The effective bias, for a particular survey, de- 
pends on the number distribution of the clusters and their 
masses, which- in turn- depends on the cosmological pa- 
rameters (mainly erg) and the cluster structure. Hence, 
constraints on dg and cluster structure from using the 
power spectrum are very degenerate. Finally, we should 
also keep in mind that cjfkp depends not only on vol- 
ume but also on number density. Hence, the area and the 
depth of the surveys affect the uncertainty estimation in 
a non-trivial manner. Constraints from adding the power 
spectrum to the cluster redshift distribution are shown by 
dot-dashed lines in Figures 2 and 3. The fact that param- 
eter degeneracies from the redshift distribution and Pd are 
quite different is evident from the rotation of the elliptical 
joint parameter constraints when one adds Pd constraints 
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to that from dn/dz. 

4.3. Including Limited Mass Follow-up 

Yet another way to reduce the cost of forcing the sur- 
vey to self-calibrate (i.e. solve for cluster structure and 
its evolution in addition to cosmology), is to add infor- 
mation from a mass follow-up of 100 clusters between the 
redshifts of 0.3 and 1.2, with 10 clusters followed-up at 
redshift intervals of 0.1. At each redshift, the 10 clus- 
ters are equally divided in mass between IO^^Ji-^Mq and 
lO^'''ft-~^M0. This is equivalent to following up less than 
1% of the cluster sample in each survey, and we envi- 
sion this followup as mostly being available through deep, 
archival X-ray observations that are being accumulated by 
XMM-Newton and Chandra even now. We see that follow- 
up significantly reduces the uncertainties for all the three 
surveys. The constraints for w goes down to ~ 0.06,0.115 
and 0.077 for SPT, Planck and DUET respectively. There 
are also major improvements in Qm (by a factor of 1.3-2.6) 
and CTg (by a factor of 2-4.5) uncertainties. Essentially, 
follow-up helps is tightening the constraints on the un- 
known cluster structure and its evolution; this translates 
into more accurate mass estimates for the clusters, and 
that reduces the uncertainties on the cosmological param- 
eters. Because SPT probes to higher redshift, follow-up is 
most valuable for this survey which contains the most high 
redshift clusters, whose masses are most uncertain without 
the follow-up information. Short-dashed lines in Figures 
2 and 3 show joint constraints when one added informa- 
tion from a follow-up of 100 clusters to that from dN/dz. 
MM02 demonstrated that although there is some further 
improvement in constraints as one follows up a larger frac- 
tion of the clusters, the difference between following up 
10% and 100% of the cluster population is minimal. We 
remind the reader that follow-up does not explicitly de- 
pend on as, n, h and fls and hence cannot constrain these 
parameters by itself and contains no cross-correlation of 
other cluster structure or cosmological parameters with 
these four parameters. 

4.4. Combining dn/dz, Pd and Follow-up 

It is interesting at this point to compare and combine pa- 
rameter constraints obtained by the two options at hand; 
follow-up mass measurements versus addition of P^i con- 
straints. Other than the Planck survey, having mass 
follow-up helps in tightening constraints on w more than 
adding power spectrum information. Planck, which sur- 
veys the maximum solid angle, delivers power spectrum 
constraints that are slightly more powerful than follow- 
up. However, for and erg, the addition of Pd con- 
straints fares better than addition of follow-up constraints 
(for example, for SPT survey addition of Pd improves the 
constraint of SIm by a factor of 2 compared to ^ 1.3 
when one adds the follow-up Fisher matrix). 

As a simple extension, one can add both power spectrum 
and follow-up constraints to that from dN/dz. This, then, 
gives us the best possible constraints on all the cosmolog- 
ical parameters. One ends up having 3.76%, 4.04% and 
6.25% constraints on w from SPT, Planck and DUET re- 
spectively. This can be compared to the overly optimistic 
3.71%, 5.05% and 3.97% constraints one obtains in the 
'Only Cosmology' case. For the DUET and Planck sur- 



veys, which are more sensitive to the massive clusters com- 
pared to SPT, the final constraints on Q.m and cts are even 
better than the 'Only Cosmology' case. Thus, one not only 
regains any cosmological sensitivity lost due to uncertain 
knowledge of cluster structure and evolution, but infact 
does marginally better. The best case scenario for pa- 
rameter constraints estimated by adding constraints from 
dN/dz, power spectrum and follow-up are shown by solid 
lines in Figures 2 and 3. In Table 1, numbers in paren- 
theses show the result of adding CMB constraints to that 
from cluster surveys. In adding constraints from CMB ob- 
servations to those from cluster surveys, one should add 
the full covariance matrix in order to get uncertainties on 
individual parameters as well as parameter degeneracies. 
Our estimate of the uncertainties in parentheses (in Table 
1) uses the CMB priors on the individual parameters and 
does not include information on their cross-correlations; 
thus, these numbers should be viewed as rough estimates. 
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Fig. 3. — Constraints on erg and Qm for an SZE SPT survey (top), 
SZE Planck survey (middle) and an X-ray DUET survey (bottom). 
Contours denote joint la constraints in four scenarios: constraints 
from dN/dz for 'Only Cosmology' case (i-dotted) ; constraints from 
dN/dz for 'Non-Std Evol' case (ii— long— dashed); constraints from 
dN/dz -\- P^i, (iii-dot-dashed); dN/dz -\- 100 cluster_follow-up, (iv— 
short-dashed) and finally constraints from dN/dz -{- P^i + fo Uow-up 
(v-solid). A flat universe is assumed for all the cases. The different 
terms are described in the text. 

4.5. Combining CMB with Pd 

Cosmological parameter constraints can, of course, be 
derived using only the cluster power spectrum Pd; how- 
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ever, taken by itself, the cluster power spectrum con- 
straints are quite weak due to large parameter degenera- 
cies. Because the shape of the power spectrum depends on 
combinations of O a/ , h and b , the best constraints are 
on these parameters and other parameters like trg and w 
are quite weak. This trend has already been noted in re- 
lated works using correlation functions (sec Rcfrcgicr ct al. 

2002) . The cluster power spectrum, by itself, is unable to 
constrain the cluster structure and evolution parameters. 
Of the three surveys, Planck constrains the cluster pa- 
rameters a few orders of magnitude better than the other 
surveys (although still extremely weakly; for example they 
are about a factor ~ 5 weaker than those from dN/dz 
alone). One can break these degeneracies by adding CMB 
anisotropy constraints and cross-correlations. Infact, if 
one puts in WMAP CMB priors on Om, erg and w, then 
one can get constraints on w of ~ 10%, although this is 
not much better than the CMB prior on w of 11%. Note 
that in contrast, once dn/dz and Pd constraints are com- 
bined, any further addition of strong external CMB priors 
on Om, f8 and w has little effecton the constraints (see 
numbers within parentheses in Table 1). In fact, it's im- 
portant to stress that future, self-calibrating cluster sur- 
veys by themselves can give roughly a factor of 2.5 (e.g 4% 
for SPT) improvement on w constraints when compared 
to those obtained using present CMB data (Spergel et al. 

2003) . 




0.24 0.28 0.32 

Fig. 4. — Constraints on Qe and f2jvf for the three surveys; SPT 
(short— dashed lines), Planck (long-dashed lines) and DUET (solid 
lines). The thin lines are for freely varying w and the thick lines are 
for fixed w = —I. 

4.6. Solving for Curvature and the Equation of State 

We have also examined the effects of our flat universe 
prior by allowing fltot to vary. This gives us an opportu- 
nity to examine the ^Ie constraints and the Qe-^m de- 
generacy (shown in Figure 4). We look at both cases of a 
freely varying w and a fixed w = — 1. We find the best con- 
straint oh^Ie to be given by the SPT survey, and these are 
= 0.0667(0.0185) when w is free {w fixed), followed by 



Planck and then DUET. When we forego the assumption 
of a flat universe and also solve for the curvature along 
with the all the other parameters, we see only marginal 
degradation in the ^Im or erg constraints. As an exam- 
ple, for Planck, the respective constraints are 0.0102 and 
0.0211 when curvature is free. On the other hand, the w 
constraints are much affected, and they degrade by factors 
greater than 3 to cr^ = 0.1353, 0.1092 and 0.1955 for SPT, 
Planck and DUET, respectively. Note that Jig constraints 
from these surveys are better than the WMAP constraints, 
and although the w constraints have weakened, they are 
also comparable to or better than constraints the cur- 
rently available tightest constraints that come from a com- 
bination of WMAP CMB measurements together with the 
galaxy power spectrum (Spergel et al. 2003) . We note that 
the parameter degeneracies are potentially quite comple- 
mentary to both CMB anisotropy and distance measure- 
ments, and thus one can obtain tighter constraints when 
one combines these different probes of dark energy. 

5. DISCUSSION AND CONCLUSIONS 

Any attempt to precisely measure cosmological param- 
eters, especially w, from future cluster surveys will re- 
quire an understanding of the mass-observable relations 
for clusters and the evolution of these relations. This is 
true whether clusters are found using X-rays, the SZE dis- 
tortion, light from the galaxies or the weak lensing shear. 
Here we examine the effect of current uncertainties about 
cluster structural evolution for three future surveys and 
require that information from the surveys be used to cali- 
bration the structure and evolution of clusters in addition 
to solving for cosmology. We show that the process of self- 
calibration weakens cosmological constraints by a factors 
of 4 — 5, even for the deepest survey, in comparison to con- 
straints derived by assuming perfect knowledge of cluster 
structure and evolution. 

We have demonstrated that adding information from 
the redshift averaged cluster power spectrum to that from 
cluster counts would help in reducing the constraints up 
to factor ~ 4. Measurement of cluster power spectrum 
helps us get a handle on the the masses of the tracer par- 
ticles through the mass dependence of the effective bias 
beff- Moreover, the cluster power spectrum delivers an 
independent measure of the parameters important to the 
shape of the transfer function- primarily U,m, h and Ob, 
and the parameter degeneracies are also different to those 
from dN/dz. Estimation of power spectrum from a clus- 
ter survey designed to measure the redshift distribution 
of clusters would not require any additional observational 
effort and hence comes rather cheaply. However, to get 
the best power spectrum constraints, one needs to survey 
large, contiguous regions of the sky. 

We have also examined the utility of detailed mass fol- 
lowup observations of 100 clusters, distributed in a red- 
shift range of 0.3 < z < 1.2. This foUowup information 
is being enabled by deep Chandra and XMM-Newton ob- 
servations today, and this information will help determine 
the parameters of cluster mass-observable relations and 
their evolution. 

The best constrains on the nature of dark energy and 
other cosmological and cluster parameters are obtained 
when the cluster redshift distribution is combined with 
the cluster power spectrum and the mass follow-up. One 
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can achieve ^ 4% constraint on w for the two SZE sur- 
veys. This is a factor of 2.5 better than the best currently 
available constraints (Spergel et al. 2003). We empha- 
size that this constraint is from a self-calibrating cluster 
survey, where the cluster mass-obscrvablc relation and its 
evolution are pulled directly from the survey and foUowup. 
Recently, Hu (2003) has shown that by including the shape 
of the mass function within survey redshift bins one can 
essentially solve for the cluster mass-observable relation as 
a function of redshift without including any information 
from the cluster power spectrum or follow-up mass mea- 
surements. Clearly, the best approach will include all the 
available information- the mass function, redshift distri- 
bution, cluster power spectrum and detailed mass mea- 
surements. 

In constructing F?-, we have subdivided each cluster 
sample into three redshift bins such that each bin has 
equal number of clusters. It is important to have thick 
redshift bins so as to avoid cross-correlation between ad- 
jacent bins. The final constraint would depend to a small 
degree on the binning of the cluster sample. For exam- 
ple, for SPT, if we subdivide the clusters in redshift bins 
z = - 0. 3, 0. 3 -0.7, 0.7 -3.0, then cr^ bom dN / dz + Pa is 
0.1577 and very close to 0.1543 which was obtained when 
the clusters were divided to have equal number in each bin. 
Incidentally, if one has a single bin from z = — 3.0 rather 
than three redshift bins, one gets au, = 0.1568. However, 
depending on the surveys, optimal redshift binning and 
the use of the full covariance matrix are important issues 
to consider in subsequent works. 

To conclude, in the present paper we have demonstrated 



the self-calibrating nature of cluster surveys and shown 
that the cluster survey is a powerful technique for studying 
the nature of the dark energy- as well as any other char- 
acteristic of the universe that affects the expansion his- 
tory, growth of density perturbations or the nature of the 
transfer function. Constraints from cluster surveys tend 
to exhibit parameter degeneracies that are quite different 
from other techniques such as CMB anisotropy studies or 
distance measurements. By pursuing all these complemen- 
tary means of studying the nature of dark energy and dark 
matter, we will produce independent constraints that will 
serve as cross checks of the different systematics underly- 
ing each technique. In addition, because the cluster survey 
is generating cosmological constraints from different phys- 
ical properties of the universe than other techniques (i.e. 
surveys probe the growth of density perturbations whereas 
distance measurements only probe the expansion history), 
the pursuit of these different techniques makes it possible 
to precisely test the current framework for understanding 
our evolving universe. We look forward to the results of 
these precise, new tests with great anticipation. 
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